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PREFACE 

Summarized  herein  are  important  findings  obtained  from  unpublished  model 
studies  conducted  by  the  Coastal  Engineering  Research  Center  (CERC)  for  the  US 
Army  Engineer  District,  Detroit  (NCE),  and  the  US  Army  Engineer  District, 
Jacksonville  (SAJ). 

This  report  was  prepared  at  the  US  Army  Engineer  Waterways  Experiment 
Station  (WES).  The  data  were  compiled  and  analyzed  by  Mr.  John  P.  Ahrens, 
Oceanographer,  and  Ms.  Martha  S.  Heimbaugh,  Civil  Engineer,  both  of  CERC.  The 
data  were  collected  by  Messrs.  Martin  Titus  and  Louis  Meyerly  and  Ms.  Karen 
Zirkel,  Civil  Engineering  Technicians,  CERC. 

General  supervision  was  provided  by  Dr.  James  R.  Houston,  Chief,  CERC, 

Mr.  Charles  C.  Calhoun,  Jr.,  Assistant  Chief,  CERC,  Dr.  Charles  L.  Vincent, 
Program  Manager,  CERC,  and  Mr.  C.  E.  Chatham,  Jr.,  Chief,  Wave  Dynamics 
Division,  CERC. 

COL  Dwayne  G.  Lee,  CE,  was  Commander  and  Director  of  WES  during  the 
preparation  and  publication  of  this  report.  Dr.  Robert  W.  Whalin  was 
Technical  Director. 
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APPROXIMATE  UPPER  LIMIT  OF  IRREGULAR  WAVE  RUNUP  ON  RIPRAP 

PART  I :  INTRODUCTION 

1.  In  many  locations  riprap  is  the  preferred  type  of  shore  protection 
against  wave  attack.  The  two  reasons  for  this  are  the  low  cost  and  nigh  dura¬ 
bility  of  stone,  and  the  effectiveness  of  randomly  placed  stone,  because  of 
its  roughness  and  porosity,  in  dissipating  wave  energy  and  attenuating  runup. 
Because  of  these  reasons,  riprap  has  been  the  most  studied  type  of  revetment, 
and  its  performance  is  well  documented. 

2.  Runup  is  one  of  the  most  important  factors  affecting  the  design  of 
revetments  exposed  to  wave  action.  Generally,  riprap  revetments  are  designed 
so  that  little  or  no  runup  exceeds  the  top  of  the  protection.  Because  of  the 
inherent  complexity  of  natural  wave  trains  and  the  interaction  of  incident 
waves  and  the  return  flow  of  previous  runup  on  a  rough,  porous  slope,  it  is 
difficult  to  predict  the  upper  limit  of  wave  uprush  on  riprap.  This  report 
summarizes  the  most  important  results  from  two  unpublished  studies,  and  pre¬ 
sents  formulas  to  calculate  the  approximate  limit  of  wave  runup.  Both  studies 
included  laboratory  tests  of  riprap  exposed  to  irregular  wave  action.  The 
formulas  can  be  used  to  compute  the  elevation  to  which  protection  needs  to  be 
extended  to  prevent  exceedance  by  runup  or  to  estimate  the  potential  severity 
of  wave  overtopping. 
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PART  II:  SOURCES  OF  DATA,  TEST  SETUPS,  AND 
TEST  CONDITIONS 


3.  The  sources  of  data  for  this  report  came  from  model  studies  conducted 
primarily  for  two  US  Army  Corps  of  Engineer  Districts,  Detroit  (NCE)  and 
Jacksonville  (SAJ). 

4.  Model  studies  conducted  for  NCE  were  to  investigate  wave  runup  on 
riprap-protected  dredge  disposal  dikes  in  the  Great  Lakes.  The  scope  of  this 
study  was  expanded  to  include  an  unusually  wide  range  of  water  depths  at  the 
toe  of  the  structure  ds  ,  zero-moment  wave  heights  Hmo  ,  and  period  of  peak 
energy  density  of  the  incident  wave  spectrum  Tp  .  By  expanding  the  scope  of 
this  study  beyond  the  immediate  problems  occurring  on  the  Great  Lakes,  the 
opportunity  to  develop  a  general  wave-runup  prediction  method  was  provided.  A 
summary  of  test  conditions  for  both  the  NCE  and  SAJ  studies  is  given  in 
Table  1,  and  data  collected  on  both  studies  are  tabulated  in  Appendix  A. 


Table  1 


Summarv  of  Test  Conditions 


1  on  2  11.9-38.5  4.9-17.5  1.02-4.74  189 

1  on  3  19.0-23.8  3.0-10.5  1.39-1.46  63-3-67.0 

1  on  4  19.0-23.8  3.2-10.3  1.39-1.49  56.9 


5.  A  1-  on  2-  (1  vertical :2  horizontal)  structural  slope  was  used  in  the 
NCE  study.  Plywood  roughened  with  glued  on  pea  gravel  was  used  as  the  sup¬ 
porting  slope  and  simulated  the  impermeable  core  of  the  dike.  This  slope  was 
covered  with  a  filter  layer  of  Sioux  Quartzite  5.5  cm  thick.  The  range  of 
weight  for  this  stone  was  from  6  to  41  g  with  a  median  weight  of  18  g.  Riprap 

armor  stone  was  placed  by  hand  on  top  of  the  filter  layer.  Armor  stones  were 

composed  of  Kimmswick  Limestone  with  a  range  of  weight  from  144  to  233  g  with 
a  median  weight  of  189  g.  The  armor  layer  was  10  cm  thick.  Figure  1  shows  a 

profile  view  of  the  model  structure,  and  Figure  2  shows  a  plan  view  of  the 


test  setup. 
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Figure  1.  Profile  view  of  structure  tested  during  NCE  study 
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Figure  2.  Plan  view  of  test  setup  for  NCE  study 


6.  The  model  was  constructed  at  a  nominal  1 : 1 6  (model rprototype)  undis¬ 
torted  Froude  scale.  The  influence  of  scale  effects  at  this  large  a  scale  was 
considered  to  be  small.  To  further  reduce  the  possibility  of  scale  effects, 
the  stone  used  in  the  filter  layer  was  double  the  size  required  for  geometri¬ 
cal  similitude.  Use  of  somewhat  larger  filter  stones  helps  establish  the 
proper  flow-regime  in  the  model  filter  layer  when  the  revetment  is  exposed  to 
wave  action;  see  Broderick  and  Ahrens  (1982).  Large-sized  filter  stone  and  a 
1:16  scale  were  used  in  both  the  NCE  and  the  SAJ  studies  to  minimize  the  in¬ 
fluence  of  scale  effects. 

7.  Tests  for  this  study  and  the  SAJ  study  were  conducted  in  a  6l-cm-wide 
channel  within  the  Coastal  Engineering  Research  Center's  (CERC's)  1.2-  by  4.6- 
by  42.7-m  wave  tank.  Wave  conditions  were  measured  offshore  by  using  three 
parallel  wire-resistance  wave  gages.  Incident  and  reflected  wave  spectra  were 
resolved  using  the  method  of  Goda  and  Suzuki  (1976).  Figure  2  shows  a  plan 
view  of  the  wave  tank  setup  for  this  study.  Details  relating  to  spectral  wave 
generation  and  the  analysis  of  wave  conditions  in  this  wave  tank  are  given  by 
Seelig  (1980). 

8.  Maximum  wave  runup  elevations  were  obtained  by  visual  observations 
made  by  an  experienced  observer,  and  quantified  by  using  a  point  gage.  The 
observer  stood  immediately  adjacent  to  the  structure  in  a  wave  absorber  chan¬ 
nel  as  shown  in  Figure  2.  The  duration  of  the  runup  observation  was  256  sec, 
corresponding  to  the  data  acquisition  system's  sampling  interval  for  the  wave 
gages  to  obtain  the  wave  information.  The  observer  tried  to  measure  the  ex¬ 
treme  excursion  of  "green"  water  near  the  middle  of  the  structure.  Observa¬ 
tions  were  not  intended  to  measure  the  upper  limit  of  spray  or  splash.  Prior 
to  using  visual  observations,  some  effort  had  been  expended  in  trying  to  use 
various  types  of  continuous  wave  gages  positioned  Just  above  the  armor  sur¬ 
face,  but  runup  elevations  that  were  measured  by  the  wave  gages  proved  to  be 
unreliable.  After  some  initial  observations  and  discussion,  two  experienced 
observers  could  obtain  maximum  runup  elevations  to  within  about  a  difference 
of  3  percent  or  less  of  each  other.  Additional  information  about  the  NCE 
study  is  given  in  Ahrens  and  Seelig  (1980). 

9.  The  SAJ  study  was  conducted  to  investigate  the  stability  of  and  wave 
runup  on  riprap  to  be  used  to  protect  Herbert  Hoover  Dike  on  Lake  Okeechobee, 
Florida.  Two  structural  slopes  were  tested  during  this  study,  1  on  4  and  1  on 
3.  Figure  3  shows  a  profile  view  of  the  1-on-4  slope  tested.  Figure  4  shows 


the  cut  and  fill  strategy  used  to  construct  a  1-on-3,  riprap-protected  slope 
on  the  embankment.  Figure  5  shows  a  profile  view  of  the  1-on-3  slope  tested 
and  the  location  of  the  wave  gages.  Figure  6  shows  a  plan  view  of  the  test 
setup, 

10.  Since  the  armor  stone  planned  to  protect  the  dike  was  marine  lime¬ 
stone  to  be  quarried  in  Florida,  this  type  of  stone  was  used  in  the  model 
tests.  This  stone  has  a  density  of  2.55  g/cm^.  The  armor  stone  had  a  median 
weight  which  ranged  from  about  57  to  67  g  during  the  course  of  the  study  (see 
Table  1).  Filter  stone  had  a  median  weight  of  about  12  g  and  a  layer  thick¬ 
ness  of  2.5  cm.  Additional  details  relating  to  test  procedures  and  setup  are 
given  in  Ahrens  and  Zirkle  (1982). 
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Profile  view  of  1-on-3,  riprap-protected  slope 
and  offshore  wave  gages 
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PART  III:  ANALYSIS  OF  DATA  AND  DEVELOPMENT 
OF  RUNUP  FORMULAS 


11.  The  biggest  difficulty  with  analyzing  the  data  from  the  NCE  study 
was  making  accurate  estimates  of  the  zero-moment  wave  height  H^  at  the  toe 
of  the  structure.  In  the  NCE  study  the  wave  heights  were  measured  offshore  in 
a  water  depth  25  cm  greater  than  at  the  toe  of  the  structure.  Due  to  shoaling 
and  breaking,  a  wide  range  of  offshore  wave  conditions  can  yield  the  same 
zero-moment  wave  height  in  shallow  water.  Therefore,  the  offshore  wave  height 
is  not  as  useful  as  the  wave  height  at  the  toe  of  the  riprap  structure.  The 
wave  conditions  near  the  structure  correlate  well  with  the  runup  and  often  can 
be  estimated  accurately  by  depth-limited  considerations.  Originally  in  the 
NCE  study  the  wave  heights  at  the  toe  were  estimated  by  using  the  method  of 
Goda  (1975)  which  accounts  for  shoaling  and  breaking  of  irregular  waves. 
However,  after  scrutinizing  the  information  generated  by  Goda's  model,  it  was 
observed  that  for  some  situations  the  method  yielded  values  of  H^/dg 
greater  than  0.8  which  is  higher  than  has  been  observed  in  any  of  CERC's  wave 
tank  calibration  tests.  Because  of  this  limitation,  it  was  decided  to  try  and 
develop  another  method  to  estimate  Hmo  at  the  toe  of  the  structure. 

12.  Several  methods  were  tested  to  account  for  the  wave  shoaling  and 
breaking  between  the  offshore  gages  and  the  toe  of  the  structure.  The  method 
that  worked  best  was  a  hybrid  method  which  combined  linear-wave  shoaling  with 
the  relation  given  by  Hughes  ( 1 984 )  as 


H 

mo 

H 

mo 

k,r\ 

i 

ksn 

where 

L  =  Airy  wave  length  calculated  at  those  depths  for  the  period  of 
peak  energy  density  Tp 

I  and  0  =  inshore  and  offshore  water  depths,  respectively 

From  wave-tank  calibration  tests  it  has  been  found  that  the  approximate 
limiting  value  for  the  zero-moment  wave  height  is  given  by 

(V)  =  0.10  tanh/^\ 


(2) 
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where  ds  is  the  water  depth  at  or  near  the  structure  toe.  The  procedure 
used  to  calculate  the  zero-moment  wave  height  at  the  toe  of  the  structure  was 
to  calculate  the  value  by  using  both  linear  shoaling  and  Equation  1  and  then 
taking  the  average  of  the  two  estimates.  If  the  average  exceeded  the  maximum 
value  suggested  by  Equation  2,  then  that  limiting  value  was  used. 

13-  The  ability  of  the  above  procedure  to  estimate  in  shallow 

water  is  demonstrated  in  Figure  7  using  wave-tank  calibration  data  collected 


PERFECT  COR  RE  LA  TION 
AND  PREDICTION  LINE  - 


OBSERVED.  H  CM 


Figure  7.  Predicted  versus  observed  for  tank  calibration  data 

prior  to  a  study  of  wave  overtopping  of  a  seawall  (Ahrens,  Heimbaugh,  and 
Davidson  1986).  This  calibration  data  included  a  wide  range  of  wave  periods 
and  an  extensive  amount  of  wave  shoaling  and  breaking  for  many  conditions 
between  the  offshore  wave  gages  and  an  inshore  gage  located  in  front  of  a  wave 
absorber  beach.  For  the  calibration  data  shown  in  Figure  7,  the  offshore 
water  depth  ranged  from  61.9  to  66.2  cm;  the  inshore  water  depth  ranged  from 
22.9  to  27.2  cm;  the  offshore  ranged  from  1.6  to  21.5  cm;  the  inshore 

Hjjjq  ranged  from  1.5  to  16.4  cm;  and  the  period  of  peak  energy  density  ranged 
from  1.75  to  3.00  sec.  The  hybrid  method  given  above  appears  to  work  well  for 
CERC  calibration  data  because  linear  shoaling  tends  to  overestimate  inshore 


Hjjjq  while  Equation  1  tends  to  underestimate  inshore  ,  and  Equation  2 

provides  a  logical  limiting  value  on  .  Based  on  the  success  of  the 

hybrid  model  in  predicting  known  data,  the  model  was  applied  to  the  NCE  study 
to  estimate  at  the  toe  of  the  structure.  Subsequent  analysis  of  pre¬ 

dicted  and  observed  maximum  runup  elevations  suggests  that  the  hybrid  method 
makes  good  estimates  of  in  shallow  water. 

14.  The  wide  range  of  water  depths  tested  in  the  NCE  study  had  been 
included  partly  to  investigate  the  influence  of  water  depth  on  wave  runup. 

This  concern  is  strongly  reflected  in  the  discussion  of  wave  runup  in  the 
Shore  Protection  Manual  (1984).  From  previous  studies  it  was  known  that  runup 
would  be  strongly  influenced  by  the  surf  condition  on  the  structures  (Ahrens 
and  McCartney  1975),  but  it  also  seemed  logical  that  the  maximum  runup  would 
be  dependent  on  the  shape  of  the  wave-height  distribution  and  nonlinear 
effects.  The  last  two  influences  would  be  very  dependent  on  the  water  depth 
at  the  toe  of  the  structure  and  the  wave  periods.  To  investigate  the  influ¬ 
ence  of  surf  characteristics  on  runup,  the  surf  parameter  for  irregular  waves 
E  is  defined  as 


tan  6 
,  \  1/2 


where 

tan  6  =  tangent  of  the  angle  8  between  the  structure  slope  and 
the  horizontal 

Lo  =  eTp  / 2-n  =  the  deep-water  wave  length 
g  =  the  acceleration  of  gravity 

When  a  runup  model  was  formulated  using  the  surf  parameter,  it  was  found  to 
contain  3ome  systematic  errors  which  could  be  related  to  the  relative  wave 
height  HjjjQ/dg  .  However,  when  a  surf  parameter  was  defined  by  using  the 
local  wave  length,  a  model  could  be  formulated  which  did  not  include 
systematic  errors  related  to  the  relative  wave  height.  The  modified  surf 
parameter  is  defined 

r  =  9 

X  ,1/2 


£ 


-/  v  v 
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where 


Lp  =  the  Airy  wave  length  calculated  by  using  the  water  depth  at  the  toe 
of  the  structure 

dg  =  the  period  of  peak  energy  density,  Tp 


Runup  is  computed  using  the  formula 

D 

max  aS 

H  =  1.0  +  bS 
mo 


(5) 


where 


Rmax  =  elevation  of  maximum  wave  runup 

S  =  surf  parameter  defined  by  either  Equation  3  or  4  depending  on 
the  prediction  method  selected 

a  and  b  =  dimensionless  runup  coefficients  determined  by  regression 
analysis 


Equation  5  has  a  form  which  is  especially  convenient  and  logical  for  predict¬ 
ing  wave  runup  on  rough  porous  slopes  as  shown  by  Ahrens  and  McCartney  (1975), 
and,  subsequently,  by  Seelig  (1980)  and  US  Army  Corps  of  Engineers  (1985).  For 
the  NCE  data,  using  the  modified  surf  parameter  defined  by  Equation  4,  the 
coefficients  in  Equation  5  were  found  to  be  a  =  1.062  and  b  =  0.153  . 

15.  By  using  the  coefficients  given  above,  Equation  5  does  a  good  Job  of 
predicting  Rmax  for  both  the  NCE  and  SAJ  studies.  Figure  8  shows  the 
predicted  values  of  Rmax  versus  the  observed  values  of  Rm  by  using  dif¬ 
ferent  symbols  to  identify  the  two  studies.  The  good  fit  to  the  NCE  data  is 
gratifying  considering  the  problem  related  to  estimating  the  Hmo  at  the  toe 
of  the  structure.  The  good  fit  to  the  SAJ  data  is  somewhat  surprising  consid¬ 
ering  that  the  thickness  of  the  armor  layer  for  the  SAJ  test  was  considerably 
thinner  than  the  armor  layer  used  in  the  NCE  tests.  In  addition,  the  struc¬ 
tural  slopes  tested  for  SAJ  were  1  on  3  and  1  on  4  compared  with  a  slope  of 

1  on  2  for  the  NCE  tests.  These  findings  indicate  that  the  maximum  runup  may 
not  be  too  sensitive  to  the  armor-layer  thickness,  and  that  the  surf  parameter 
properly  accounts  for  differences  in  the  structural  slopes.  It  should  also  be 
recalled  that  the  runup  coefficients  were  obtained  from  the  NCE  study  so  that 
the  SAJ  data  provide  a  rather  severe  test  for  the  runup  model's  predictive 
ability. 

16.  By  lumping  the  data  from  the  two  studies  together,  somewhat  better 
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Figure  8.  predicted,  using  coefficients  from  NCE  data,  versus 

Rmax  observed,  both  NCE  and  SAJ  studies 


a  and  b  runup  coefficients  can  be  determined.  By  using  regression  analysis 
on  the  combined  data  set,  the  improved  coefficients  are  a  =  1.154  and 
b  =  0.202  .  A  new  scatter  plot  comparing  predicted  and  observed  values  was 
prepared  with  the  above  coefficients,  and  Equation  5  was  used  to  calculate  the 
predicted  values  of  maximum  runup.  The  new  scatter  plot  (see  Figure  9)  shows 
that  the  change  in  the  runup  coefficients  caused  very  little  change  over  the 
scatter  plot  shown  in  Figure  8.  Even  though  there  was  little  change  in  the 
scatter  plot,  the  limiting  value  for  R/H^  dropped  from  6.9  to  5.7;  the 
limiting  value  for  Equation  5  is  given  by  the  ratio  of  a  to  b  . 

17.  To  investigate  systematic  error  in  predicting  the  maximum  runup  and 
to  identify  possible  ways  to  improve  the  prediction  method  based  on  the  modi¬ 
fied  surf  parameter  c  ,  a  series  of  error  plots  was  made.  In  these  plots, 
the  percent  error  %E  in  predicting  the  maximum  runup  is  defined  as 
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Figure  9. 


Rmax  predicted,  using  coefficients  from  both  NCE  and  SAJ  data, 
versus  observed  for  NCE  and  SAJ  data 


(Rmax)p  ~  (Rmax)t 


where  the  subscripts  p  and  o  indicate  predicted  and  observed,  respec¬ 
tively.  The  percent  error  is  plotted  versus  ,  J-  ,  ds/Lp  *  Hmc/^p  * 

H^/dg  ,  r(bar)/d^0  ,  and  II,  and  cot  6  in  Figure  10a,  b,  c,  d,  e,  f,  g, 
and  h,  respectively,  where  r  is  the  average  armor-layer  thickness  and  II 
is  Goda's  (1983)  nonlinear  parameter  defined  for  irregular  waves 


tan  h- 


/2itd 
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The  larger  the  value  of  II  ,  the  more  nonlinear  the  waves  with 
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—  ,  for  deep-water  conditions 
P 


II-*-  ■  °°— P  —  ,  Ursell's  parameter  for  shallow-water  wave  conditions 

K)3 

Figure  10  shows  little  or  no  systematic  error  for  the  prediction  method  based 
on  Lp  .  Figure  10  also  shows  that  the  percent  error  ranges  from  -33  to 
+44  percent,  but  that  for  most  tests  the  error  is  within  about  ±10  percent. 

18.  Approximately  25  percent  of  the  tests  had  a  percent  error  greater 
than  ±10  percent.  Because  of  this,  it  may  be  useful  in  some  critical  or  life- 
threatening  situations  to  use  a  value  of  Ru,ax  greater  than  the  expected 
value  produced  by  Equation  5  when  using  the  recommended  coefficients.  Fig¬ 
ure  11  shows  how  the  percent  error  which  has  been  normalized  by  the  standard 
deviation  of  the  data  set  a  seems  to  have  the  shape  of  a  normal  distribu¬ 
tion.  To  test  this  hypothesis,  namely,  that  the  percent  error  has  a  normal 
distribution,  a  Kolmogorov-Smirnov  (K&S)  test  was  performed.  This  test  is 
used  to  determine  whether  or  not  the  data  deviate  a  statistically  significant 
amount  from  the  assumed  normal  distribution  model  (Cornell  and  Benjamin  1970). 

19.  The  K&S  test  indicates  that  the  normal  distribution  for  error  should 
be  accepted  at  the  20-percent  significance  level.  A  20-percent  level  is  a 
more  severe  criterion  than  a  iO-percent  level  as  it  indicates  there  is  a 
20-percent  chance  of  rejecting  a  model  which  is  in  fact  true,  a  Type  I  error. 
The  20-percent  significance  level  is  the  most  severe  criterion  commonly  tabu¬ 
lated  for  the  K&S  tests.  Recognizing  that  errors  have  a  normal  distribution 
provides  an  easy  way  to  give  more  conservative  estimates  of  Rmax  than  is 
provided  by  a  regression  equation.  Generally,  about  half  the  errors  are  above 
the  regression  curve  and  about  half  are  below,  so  the  curve  represents  a  50_ 
percent  exceedance  level.  In  Figure  12  a  more  conservative  trend  is  shown 
above  the  regression  curve.  The  conservative  curve  was  generated  by  increas¬ 
ing  the  runup  regression  coefficient  a  by  two  standard  deviations  of  the 
percent  error,  i.e., 

ac  (conservative  a)  =  a  (1  +  2o)  =  1,143  (1.0  +  2  x  0.1286)  =  1.437 
The  value  of  the  runup  coefficient  b  remains  the  same,  i.e.,  b  =  0.202  . 
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It  can  be  seen  in  Figure  12  that  the  new  conservative  runup  curve  provides  a 
envelope  for  the  observed  data.  The  conservative  curve  would  be  expected  to 
exceed  about  97.7  percent  of  the  data.  An  exceedance  level  of  97.7  percent  is 
obtained  from  a  standard  normal  distribution  table  for  a  value  two  standard 
deviations  greater  than  the  mean.  Figure  12  helps  confirm  the  method  of 
choosing  an  envelope  curve  by  showing  only  one  observed  value  above  the  con¬ 
servative  curve.  This  is  approximately  what  would  be  expected  for  a  normal 
distribution  with  a  sample  size  of  69.  Other  curves  used  to  predict  maximum 
runup  could  be  constructed  which  would  be  more  or  less  conservative  than  the 
example  Just  provided.  The  degree  of  conservatism  would  be  evaluated  on  the 
basis  of  the  risk  posed  by  waves  overtopping  the  revetment. 

20.  Since  the  standard  surf  parameter  has  been  frequently  used  to  pre¬ 
dict  wave  runup,  it  is  useful  to  provide  a  prediction  formula  based  on  that 
method  to  allow  comparison  to  earlier  studies.  Using  Equation  3  to  define  the 
surf  parameter  in  Equation  5,  the  runup  coefficients  were  determined  for  the 
combined  NCE  and  SAJ  data  sets  as  a  =  1.022  and  b  =  0.247  .  Figure  13 
shows  the  predicted  and  observed  values  of  versus  the  standard  surf 

parameter.  It  can  be  seen  that  the  predicted  values  follow  the  trend  of  the 
observed  data  very  well.  Using  the  same  method  to  evaluate  errors  as  was  used 
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Figure  13.  Comparison  of  the  prediction  of  Rmax  using  f; 
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for  the  Equation  4  model,  it  was  found  that  there  were  no  systematic  errors 
associated  with  the  model  which  used  the  standard  surf  parameter.  Figure  14 
shows  that  the  standarized  percent  errors  for  this  model  also  seem  to  have  a 
normal  distribution.  Performing  the  K&S  test  once  again  showed  that  these 
data  were  also  normal  at  the  20-percent  significance  level  and  thus,  could  be 
assumed  to  have  a  normal  distribution.  Figure  15  shows  the  more  conservative 
curve  which  could  be  expected  to  envelop  97.7  percent  of  the  data  and  repre¬ 
sents  an  increase  of  two  standard  deviations  over  the  expected  mean  curve. 

The  coefficients  for  this  conservative  curve  are  a  =  1.285  and  b  =  0.247  . 
Once  again,  this  curve  is  only  one  of  many  more  conservative  curves  that  could 
be  constructed  depending  upon  the  design  situation.  Using  the  runup  coeffi¬ 
cients  with  the  standard  surf  parameter  would  be  an  easy  way  to  estimate 
using  a  small  calculator.  The  more  accurate  model  would  require  the  calcula¬ 
tion  of  Lp  for  use  in  the  modified  surf  parameter  which  would  be  more 
difficult  than  calculating  the  deep-water  wave  length  for  the  standard  surf 
parameter . 
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Figure  14.  Normal  distribution  curve  for  £  model  data 
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PART  IV:  SUMMARY 


21.  All  the  equations  presented  within  this  report  were  developed  from 
two  unpublished  laboratory  studies.  These  equations  provide  an  easy  way  to 
calculate  the  Rmax  of  irregular  waves  on  riprap-protected  embankments. 
Table  2  summarizes  the  important  information  about  two  runup  equations  which 

Table  2 

Summarized  Information  for  Maximum  Runup  Formula.  Equation  5 


Formula 

Category 

Wave 

Length 

Used 

Surf 

Parameter 

Used  in 
Equation  5 

Runup 

Coefficients 

Variance 

Explained 

o  ,  Standard 
Deviation  of 
Percent 
Error 

Recommended 

lp 

(Equation  4) 

a  =  1.154 

b  =  0.202 

R2  =  0.843 

12.3 

Alternative 

Lo 

5 (Equation  3) 

a  =  1.022 

b  =  0.247 

R2  =  0.817 

12.9 

represent  the  most  accurate  existing  method  to  determine  the  approximate  upper 
limit  of  wave  uprush  on  a  riprap  revetment.  The  two  equations  are  presented 
as  a  recommended  method  and  an  alternate  method  to  compute  Rmax  .  The  recom¬ 
mended  method  has  little  or  no  systematic  error  such  as  might  be  associated 
with  the  influence  of  water  depth  or  nonlinear  effects  and  is  slightly  more 
accurate  than  the  alternate  method.  The  alternative  method  is  easier  to  cal¬ 
culate  and  can  serve  as  a  "rule  of  thumb"  estimate.  In  Table  2  the  runup  co¬ 
efficients  are  to  be  used  in  the  general  runup  equation  (Equation  5)  by  using 
either  the  standard  or  modified  surf  parameter  as  noted.  A  method  was  devel¬ 
oped  which  provides  a  reasonable  way  to  make  the  predicted  values  of  R,,,^ 
more  conservative.  It  was  found  that  the  errors  in  predicting  Rmax  have  a 
normal  distribution,  and  this  fact  was  used  to  adjust  the  runup  coefficient 
a  so  that  any  predetermined  exceedance  level  for  Rmax  could  be  achieved. 

For  example,  by  increasing  the  coefficient  a  by  two  standard  deviations  of 
the  percent  error  gives  R^,^  predictions  which  would  be  expected  to  exceed 
97.7  percent  of  the  observed  values  of  Rmax  .  This  technique  produces  a 
logical  envelope  for  the  data.  Table  2  lists  the  standard  deviation  o  of 


the  percent  error  which  was  used  in  this  method  and  the  correlation  squared 
which  is  the  variance  explained  by  the  regression  equation  used  to  predict 

^max^^mo  ' 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 

22.  Visual  observations  of  the  maximum  runup  of  irregular  waves  on  rip¬ 
rap  give  reliable  values  of  the  maximum  elevation  of  wave  uprush.  For  the  two 
studies  considered  the  time  interval  of  observation  was  256  sec.  It  is  diffi¬ 
cult  for  an  observer  to  maintain  adequate  concentration  on  the  runup  process 
for  intervals  longer  than  256  sec.  This  interval  provides  between  100  and 
250  runup  events  at  the  wave  periods  tested.  Future  tests  will  consider  using 
photogrametric  methods  to  measure  irregular  wave  runup  on  riprap  to  increase 
the  time  interval  of  observation  and  to  obtain  the  entire  runup  distribution 
rather  than  Just  the  maximum  value. 

23.  The  runup  equations  presented  appear  to  be  the  best  available  to  es¬ 
timate  the  approximate  upper  limit  of  irregular  wave  uprush  on  riprap  revet¬ 
ments.  Further  tests  are  planned  which  should  produce  improved  methods  to 
determine  the  runup  characteristics  of  irregular  waves  on  rough  and  porous 
slopes. 
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Table  A1 

Summary  of  Jacksonville  (SAJ)  and  Detroit  District  (NCE)  Test  Conditions 


offshori  uvt  of fshore  inshore  estieite  inshore 

depth  period  H»o  depth  inshore  nvi 

Tp  d  H*o  length 


No. 

ce 

sec 

Cl 

Cl 

Cl 

Cl 

Cl 

_ 

34.911 

2.753 

5.177 

11.983 

4,943 

294.171 

11.392 

2 

34.911 

1.533 

6.797 

11.911 

6.394 

139.919 

9.392 

3 

34.911 

1.128 

7.340 

11.911 

6.494 

181.714 

9.514 

4 

34.911 

3.821 

18.218 

11.981 

7.544 

411.311 

24.413 

5 

34.911 

2.132 

9.536 

11.911 

7.334 

215.118 

14.897 

6 

34.911 

1.684 

13.171 

11.981 

7.212 

176.742 

28.413 

7 

43.5M 

3.282 

5.231 

2I.5M 

5.199 

459.122 

12.412 

9 

43. 5M 

1.842 

7.424 

21. 3M 

6.797 

132.112 

12.192 

9 

45.3* 

1.571 

7.788 

21.588 

7.252 

210.219 

13.584 

it 

45.381 

1.954 

9.158 

2I.5M 

8.644 

267. 8M 

17.793 

n 

45.5* 

3.587 

9.834 

20. 3M 

9.584 

491.389 

28.3* 

12 

45. 5* 

2.132 

10.710 

28.5M 

10.239 

278.342 

17.282 

13 

43. 311 

2.132 

13.402 

20. 3M 

12.274 

278.342 

22.498 

14 

43.3M 

1.684 

13.777 

20. 5M 

22.890 

227.174 

22.993 

13 

33.311 

3.202 

5.747 

28.381 

5.381 

338.626 

16.497 

16 

33.311 

1.141 

8.289 

28.381 

7.638 

143.883 

14.387 

17 

33.311 

1.381 

12.844 

28.311 

11.411 

243.884 

19.393 

19 

33.311 

1.940 

13.M4 

28.311 

14.339 

387.701 

33.393 

19 

33.311 

1.420 

17.778 

28.311 

15.691 

251. M2 

26.194 

21 

33.311 

2.170 

14.432 

28.511 

14.040 

347.837 

26.483 

21 

33.311 

4.741 

12.487 

28.381 

12.388 

785.372 

31.994 

22 

33.381 

3.281 

10.737 

28.381 

11.421 

324.679 

24.994 

23 

3B.299 

3.122 

3.924 

33.299 

5.733 

358.933 

18.891 

24 

38.299 

1.243 

8.155 

33.299 

7.649 

191.843 

13.412 

23 

38.299 

1.533 

18.331 

33.299 

9.779 

238.339 

17.697 

26 

31.299 

2.170 

14.971 

33.299 

14.394 

373.350 

25.794 

27 

38.299 

1.620 

18.419 

33.299 

17.484 

247.762 

31.884 

28 

38.299 

3.241 

9.763 

33.299 

9.486 

572. 8M 

27.794 

29 

58.299 

4.741 

15.573 

33.299 

15.214 

047.477 

37.718 

31 

38.299 

2.989 

11.881 

33.299 

10.747 

311.541 

27.794 

31 

43.388 

2.753 

6.239 

38.581 

6.147 

316.442 

17.612 

32 

43. 3M 

1.249 

9.224 

38.388 

8.704 

282.284 

19.212 

3J 

43. 3M 

1.515 

8.388 

38.588 

7.803 

240.949 

1J.584 

34 

63. 5M 

1.814 

9.670 

38.3M 

9.249 

323.817 

21.413 

35 

43.388 

1.584 

12.94* 

38. 3M 

12.281 

259.052 

22.898 

36 

63. 3M 

4.741 

14.381 

38. 3M 

14.031 

918.131 

33.989 

37 

63.3M 

1.620 

14.931 

30  581 

14.211 

283.441 

24.898 

38 

43.518 

2.813 

0.992 

38. 3N 

8.728 

521.574 

23.118 

39 

34.981 

1.062 

6.683 

11.911 

6.079 

116.411 

9.611 

41 

34.981 

3.468 

7.085 

11.981 

7.523 

371.134 

18.882 

3 

48.112 

1.491 

NA 

23.812 

7.544 
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9.611 

4 

48.112 

1.440 

NA 

23.112 

3.884 
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5.313 

a 

48.112 

1.428 

NA 

23.812 
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9.984 

9 

41.812 
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NA 

23.812 

3.181 

288.8M 
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ii 

44.838 

1.398 

NA 

19.838 

11.248 

177.888 

11.497 

u 

44.138 

1.428 

NA 

19.058 

9.438 

111.  4M 

11.314 

12 

44.838 

1.428 

NA 

19.038 

8.113 

111.  BM 

0.592 

13 

44.838 

1.470 

NA 

19.131 

6.811 

189. 4M 

7.211 

23 

48.812 

1.438 

NA 

23.812 

8.723 

284. 7M 

13.5* 

26 

48.112 

M68 

NA 

23.012 

6.411 

215.981 

9.792 

27 

48.112 

1.428 

NA 

23.012 

3.143 

199. 7M 

7.1* 

31 

48.812 

1.428 

M 

23.817 

8.882 

198,9* 

13.212 

31 

44.838 

1.410 

NA 

19.038 

18.312 

181. 3M 

12.7* 

32 

44.138 

1.398 

NA 

19.838 

9.639 

177. 4N 

11.611 

33 

44.838 

1.460 

NA 

19.838 

8.323 

184. 8M 

11.592 

34 

44.838 

1.410 

NA 

19.038 

6.218 

188.3M 

9.611 

33 

44.138 

1.468 

NA 

19.038 

2.991 

186.888 

5.611 

34 

44.838 

1.398 

NA 

19.031 

11.316 

177. 7M 

13.297 

37 

44.838 

1.428 

NA 

19.038 

9.544 

181.  SN 

12.112 

IB 

14.838 

1.418 

NA 

19.038 

8.192 

178. IN 

11.697 

13 

48.812 

1.468 

NA 

73.817 

8.648 

287.  IN 

13.5* 

44 

48.812 

1.441 

HA 

23.812 

6.494 

282.6M 

11. *1 

r> 

4B.I12 

1  441 

NA 

23.017 

3.124 

283.  IN 

6.196 

*6 

44.138 

1.448 

NA 
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183. 4M 

12.992 

47 

44.838 

J.460 

NA 
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(8 
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1.408 
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179.6* 
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44.138 
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NA 
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NA 
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